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Abstract. We search for possible time lags caused by quantum gravitational 
(QG) effects using gamma-ray bursts (GRBs) detected by INTEGRAL. The 
advantage of this satellite is that we have at our disposal the energy and arrival 
time of every detected single photon, which enhances the precision of the time 
resolution. We present a new method for seeking time lags in unbinned data using 
a maximum likelihood method and support our conclusions with Monte Carlo 
simulations. The analysis of the data yields a mass scale well below the Planck 
mass whose value may however increase if better statistics of GRBs were available. 
Furthermore, we disagree with previous studies in which a non-monotonic function 
of the redshift was used to perform a linear fit. 
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1. Introduction 

There is a general agreement that the classical space-time structure as described by the 
theory of general relativity will undergo drastic modifications at very small distances 
and very large energies due to quantum fluctuations. It is commonly argued that the 
relevant scales at which some new phenomena caused by quantum gravity (QG) occur 
are determined by a combination of Newton's constant G, Planck's constant h and the 
velocity of light in vacuo c, i.e. by the Planck length l P = yJhGjf? w 1.6 ■ 10~ 33 cm 
or equivalently the Planck energy Ep or the Planck mass Mp — Ep/c 2 — y/hc/G « 
1.2 ■ 10 19 GeV/c 2 . Until recently it was thought it would be almost impossible to 
detect the effects of such extremely short length scales or large energies. 

Although a full quantum theory of gravity has not yet been established it has been 
realized that some generic predictions seem to emerge from the various approaches to 
a theory of QG. Assuming that QG possesses a well defined semiclassical limit which 
is obtained for weak gravitational fields and/or low energies, E -C Ep, one can look 
for falsifiable predictions from semiclassical QG to first order in E/Ep. 

One of the most striking predictions is a distortion of the photon dispersion 
relation 



where E, p denote the photon energy and momentum, respectively, a is a model- 
dependent dimensionless parameter of order unity and c is the (standard low-energy) 
velocity of light in vacuo. The non-standard dispersion relation (fTJ) leads to an energy- 
dependent velocity of light, v = v(E), defined by the group velocity v := dE/dp: 



The corrections to the velocity of light of the form ([2]) could be interpreted as an 
explicit violation of Lorentz invariance at the Planck scale. For example, there may 
exist a preferred frame which is commonly chosen to be the frame that coincides with 
the rest frame of the cosmic microwave background radiation, implying that light 
would have a helicity-dependent velocity. In 1 a more general dispersion relation 
with Lorentz symmetry breaking terms that depend explicitly on the hclicity of the 
photon was studied and bounds on the QG scale from different astrophysical sources 
were given. However, observations of synchrotron radiation in the Crab nebula [2] 
tend to rule out a helicity-dependent velocity of light. 

Another possibility to understand a possible violation of Lorentz invariance, 
as proposed e.g. in non-critical string theory [3J 0], string theory [5] or effective 
field theory approaches 6 , is to interpret the energy-dependent velocity of light as 
v(E) = c/n(E), where n(E) is the refraction index of the non-trivial optical properties 
of the " foamy" structure of space-time caused by quantum fluctuations on short time 
and distance scales. 

A very promising approach to QG is Loop Quantum Gravity (for reviews see 
[3 El En [11]), where modifications of the type (P)-© are present in the 2+1 
dimensional theory and where it is conjectured [T51 HI] that the same will be 
true in QG in 3+1 dimensions. In this theory, the corrections are understood as 
indicating not a breaking of Lorentz invariance but rather a deformation of it. One 
assumes that the relativity of inertial frames is preserved, however, one requires that 
there be two constant scales which are observer-independent: the standard velocity of 



E 2 = p 2 c 2 + a— + 0(E^lE 2 p), 



(1) 




(2) 



Lorentz violation in INTEGRAL j-ray bursts 



3 



light c and the Planck length lp (or equivalently Ep). It has been shown that Lorentz 
invariant theories satisfying these requirements exist if the Lorentz transformations are 
treated not in the standard way, but are realized non-linearly when acting on energy 
and momentum eigenstates. Such theories are called Doubly or Deformed Special 



Relativity (DSR) [I3[T![T3[Ta[Ta[I3[2Tl[22[2^ 



Different realizations of DSR lead to a different energy dependence of the velocity of 
light. However, a common feature of all DSR models is that the velocity of light does 
not depend on helicity. 

Note, however, that in [32 it was argued that, in order to construct a quantum 
field theory that consistently incorporates DSR, it should not depend on extensive 
quantities like the total four momentum of particles but rather on intensive quantities 
like the fields' energy and momentum densities. This not only solves the soccer- 
ball problem but also changes drastically the predicted effects by many orders of 
magnitude. As shown there, the effect of the "new DSR" is about 57 orders of 
magnitude smaller than predicted by the "old DSR", thus making DSR hardly 
measurable. 

In view of the fact that there exists a large variety of approaches to QG which 
lead to an energy-dependent velocity of light of the form ((2]), it seems worthwhile to 
seek experimental tests of Equation ([2|) . It was pointed out that one powerful way to 
probe Equation @ may be provided by gamma-ray bursts (GRBs) [33J IMJ- Several 
studies have been conducted using measurements of GRBs [35j [36l [37] . GRBs are the 
most distant variable astrophysical sources of energetic photons detected by present 
experiments in the energy range from keV to GeV. 

In this paper we shall consider the following form of the velocity of light 



where we have put the constant a equal to ±1 and replaced the Planck mass Mp 
by a QG mass M to be determined or constrained by GRBs. Furthermore, we have 
neglected the higher order terms because the energy E of the photons emitted by the 
available GRBs detected by INTEGRAL are much smaller than the expected energy 
scale Mc 2 representing the QG effects. 

Light propagation from GRBs is not only determined by the velocity of light ([3]) 
but is also affected by the cosmological expansion of the universe. Present observations 
are consistent with a nearly (spatially) flat universe. In the following we shall assume, 
for simplicity, an exactly flat universe described by the ACDM model consisting of 
baryonic matter (bar), cold dark matter (cdm) and a positive cosmological constant 
A, i.e. f^tot = + with f2 m := f^bar + ficdm = 0.27. The time delay between two 
photons with an energy difference AE is then given by 



where Hq = 71 kms -1 Mpc -1 is the Hubble constant (Hq 1 — 13.77 Gyr). Assuming 
an energy difference AE = 300 keV and a redshift z — 3 we get a time lag of 
approximately At = 2 • 10 ~ 5 s for M = Mp. 

In spite of the fact that GRB signals are interesting for searching for QG effects, 
they are far from being perfect, mainly due to our lack of knowledge of the internal 
physical processes which are at the origin of the light emission. It is conceivable 
that photons of different energies are produced by different mechanisms within the 




(3) 




(4) 
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GRB, thus narrowing the energy range in which a comparison between arrival times 
is possible. 

In this paper we study possible time lags from GRB light curves detected by 
INTEGRAL. Contrary to previous studies [351 GS3 [33 SHI HO] where a binning in time 
and energy was used, we use unbinned data, i.e. we know the arrival time and energy 
of every single photon detected by the satellite. In order not to destroy this valuable 
piece of information we will not use wavelets as was done in the cases described above 
but rather use a new method. 

The paper is organized as follows. In scction[2]we describe the relevant properties 
of INTEGRAL and in section[3]our method to analyze the unbinned data. In section|4] 
we present the results of the Monte Carlo simulations. Section [5] contains the results 
using GRBs detected by INTEGRAL, and we finish with the conclusions in section [6] 

2. INTEGRAL satellite 

INTEGRAL [51] is a mission of the European Space Agency (ESA) devoted to gamma 
ray astronomy. It features a coded mask instrument ISGRI 52 . This instrument 
enables us to measure for each photon in the energy range 15 keV to 1 MeV the 
arrival time with a precision of 6 • 10 -5 s as well as the energy with a precision of 10%. 

The detector has a dead time of about 25%. This dead time is a function of the 
incoming rate and can vary during a GRB. The dead time is measured internally by 
the instrument and is given as a mean dead time over 8 seconds independently for 
6 parts of the detector. It can be corrected statistically in weighing each incoming 
photon by 1/(1 — dead time) with the corresponding time slice and detector part dead 
time. If the rate exceeds telemetry capabilities a data gap is created in wich the dead 
time is 100%. In this case it cannot be statistically corrected and we have a hole in the 
data versus time. This unfortunately happens frequently during very intense GRBs. 

The instrument also registers an important rate from the background due to 
diffuse photons from the sky, internal radioactivity of the instrument and flux from 
sources present in the field of view. This background rate varies with time but not 
perceptibly during the typical time scale of a GRB. We have two ways of predicting 
this background. Before and after the GRB the background can be measured as 
the full rate registered by the instrument. During the GRB, the pixels that are in the 
shadow of the mask for the direction of the GRB register only the background photons 
of the GRB. The illuminated pixels register this background as well as the flux from 
the GRBs. Statistically the rate from the GRB can be computed by properly weighed 
subtraction. As, most of the time, the GRBs are in the partially coded field of view, 
the number of pixels available for background measurement is bigger than the number 
of pixels seeing the source. 

The fraction of a pixel that is illuminated by the GRB (so called PIF value) can 
be calculated with the knowledge of the coordinate of the GRB and the knowledge of 
the attitude of the instrument. We are not able to determine individually if a photon 
comes from the GRB or the background, but the PIF can be used to properly weigh 
its probability to come from the GRB. For example, a light curve can be built by 
using only pixels that are fully illuminated by the source and removing the constant 
rate measured by the completely opaque pixels. 
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3. Description of the analysis method 

The majority of the GRBs seems to follow a pattern called Fast Raise and Exponential 
Decay (FRED). In order to model a GRB light curve, we parameterize it with five 
parameters and call the resulting probability distribution / = f(ti, Ef, P, B, R, D, k, h) 
(see figure [T]). We suppose that a set of measured parameters U and Ei came from 
the probability density function /. We use the method of maximum likelihood, which 
consists of finding the set of values P, B, R, D, k and h, which maximizes the joint 
probability distribution for all data, given by 

F(P, B, R, D, K,h)=Y[ f(ti, Ef,P, B, R, D, k, h) (5) 

i 

together with the constraint 

ti 

dt' f(t',Ei;P,B,R,D,K,h) = l, (6) 



where T is the likelihood function and the integral runs between to and t\ as shown in 
figure [TJ In fact, the condition ([6]) that the integral over time be equal to one reduces 
the degrees of freedom for / and F by one. For example, B can be chosen to be fixed 
by this condition, so we can think of / and T as not depending on B. However, for 
clarity we write the £?-term dependence for both functions. 

It is easier to search for the parameters that maximize lnJ 7 , as the products 
on the right hand side of Equation ([5]) is now a sum. To find these parameters, we 
use a multidimensional unconstrained nonlinear minimization where we minimize the 
function — In T . 



i. 




Figure 1. Sketch of a typical light curve of a GRB for a given energy interval. 
The curve is parameterized by five parameters: B is the background level, _R the 
duration of the rise, h the height above the background, D the decay time for 
exp(— t/D) and re describes the magnitude of the dependence on the energy of 
the distribution /, t m ax = P + re • E, where P is the time when the intensity 
reaches a maximum and E is the photon energy. The area under the curve must 
be one, so that one parameter, e.g. B , is fixed by this condition. The dashed 
line shows a distribution for another energy interval that is shifted by an amount 
of At = re ■ AE sketching the shift in time due to quantum gravitational effects. 
This shift is usually much smaller than the other parameters. 



Figure Q] shows a typical light curve of a GRB. We always choose time intervals so 
that such a sketch can be found. However, in order to avoid wrong results, we also take 
account for other possibilities when for example R > t mli ^(E) — to or t\ < t max (E). 
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The maximum shift in time due to quantum gravity is expected to be of the order of 
2 • 10 -5 s, which is smaller by a factor of three than the time resolution of INTEGRAL. 
Therefore, it is at first highly questionable whether such time differences can be 
measured, not to speak of the results gotten from unbinned data. In order to get a 
better feeling of the behavior of the likelihood, we performed Monte Carlo simulations 
with a total number of photons ranging from 500 to unrealistic 300'000. First, we 
created N events i with energy Ei distributed according to a typical GRB event. 
That is, a typical energy distribution for the photons of GRBs follows the pattern of 
the so-called Band function [iSS. given by the following equation: 



N E (E) = A 



A 



E 



100 keV 



(a - (3)E 



ex P 77~ 



E 



(a - (3)E > E, 



100 keV 



-0 



E 







exp(/3 — a), 



100 keV, 

(a - P)E < E, 

where we choose typical values for the parameters, i.e. a = — 1, 
E = 200 keV (see left panel of figure EJ). 
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Figure 2. Left panel: example of a photon distribution as a function of the 
energy. The energy ranges from 20 keV to 300 keV according to the Band function 
(0, the total photon number is 5000. Right panel: Example of a simulated GRB 
for k = -10 -5 s/keV, P = 1 s, R = 0.3 s, h = 50 s" 1 and D = 0.5 s with a total 
photon number of 5000. In order to be able to compare the fit with the GRB, 
we require that the areas under both curves be equal, so that the parameter B is 
recovered. The overlaid curve is the FRED function with fitted parameters for a 
photon of energy 0. As the parameter n is very close to this curve represents 
well the family of FRED curves of the problem. 



With this energy distribution, we created arrival times for each photon according 
to the FRED distribution /. In addition, because the time resolution of INTEGRAL is 
6.1 • 10~ 5 s, we perturbed the arrival time of each photon with a Gaussian distribution 
with a deviation of 6.1 • 10 -5 s. The Monte Carlo simulations were done with 
K = -10~ 5 s/kcV, P = 1 s, R = 0.3 s, h = 50 s^ 1 and D = 0.5 s. Remember 
that At = k ■ AE, so that a value for k of 10 -5 s/keV represents a maximum time 
delay of ~ 3- 10~ 3 s, which is well longer than the expected time delay due to quantum 
gravitational effects. 
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Figure [2] (right panel) gives an example of a simulated GRB for parameters as 
described above. The histogram shows a typical simulation of a GRB using a FRED 
distribution, while the black line shows the solution of the minimization of Equation 
©. This curve is defined by P = 0.983, k = -4.95 • 10" 5 s/keV, R = 0.276 s, h = 48 
s _1 and D = 0.485 s. Except the value k which is five times too big, the other values 
are easily recovered by the minimization of Equation [[5"]), However, the Monte Carlo 
simulations have a tendency to underestimate the parameters. As can be seen in 
table [H except the mean value of D for N = 500, all values are too low for small N. 
Note that, apart from k, R is not well estimated and has therefore a big deviation. 

Table 1. Results of 200 Monte Carlo simulations for each value of N 



N 




P [s] 


k [s/keV] 


R[s] 


h [s- 1 ] 


D [8] 


500 


I- 1 
a 


0.95 
0.10 


-4.7- 
4.0 • 


10" 5 
10" 4 


0.17 
0.21 


4.57 
1.4 


0.65 
0.31 


1000 


Ji 

17 


0.97 
0.08 


-6.5 • 
2.7- 


io- 7 

10" 4 


0.19 
0.21 


4.8 
1.11 


0.59 
0.24 


2000 


17 


0.99 
0.03 


-2.3- 
1.6 • 


10" 5 
10" 4 


0.19 
0.23 


5.1 
0.52 


0.51 
0.09 


5000 


A* 

(7 


0.99 
0.02 


-3.3 
8.6 • 


10" 6 
10" 5 


0.17 
0.25 


5.0 
0.28 


0.50 
0.03 


10000 


ft 
a 


1.0 
0.01 


-2.0 
6.8 • 


10" 5 
10" 5 


0.17 
0.25 


5.0 
0.2 


0.50 
0.02 


3 ■ 10 5 


A* 
a 


1.00 
0.002 


-1.0 ■ 
1.2 • 


10" 5 
10" 5 


0.30 
0.002 


5.00 
0.04 


0.50 
0.004 



From table Q] it should be clear that even with 3 • 10 5 photons it is not possible to 
get a trustful result for that small a value K. Recall that k = 10 -5 s/keV is about a 
factor 100 larger than the expected time lags caused by quantum gravitational effects. 
A crude way of evaluating the statistics necessary for a convincing measurement is 
to make the assumption that the FRED distribution may be approximated by a 
Gaussian distribution. This distribution is obtained by minimizing the error of N 
independent measurements, where the single parameters are p, and a. The error of 
a single measurement is given by cr/v^/V, so that if we want to reach a precision of 
At = 10~ 5 s with a burst lasting one second, we need 10 10 photons. 

A more careful analysis shows that the standard deviation for a FRED 
distribution does not behave like const/ yN. Figure [3] shows the standard deviation 
o~ K as a function of the photon number N as given in table [TJ A fit to the data points 
between N = 500 and N = 10000 shows that the standard deviation of a FRED 
distribution is given by 

a K = 0.0182 ■ TV" ' 617 , (8) 

where the exponent is smaller than the usual 1 j\J~N for a Gaussian. With this equation 
we are also able to assess the error for k when using data from GRBs measured by 
INTEGRAL. 
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Figure 3. Standard deviation a K for function of the number of photons 

N . The data points are shown in table [T] in the column re. The solid line shows 
the fit and is given by Equation J8}. 



5. Results from GRBs detected by INTEGRAL 

5.1. Determination of the parameter k 

The data provided by INTEGRAL contains for each single registered photon four 
pieces of information: the arrival time, the energy, the dead time and the PIF value 
(see section [5]). In our analysis we take only photons that have a PIF value larger 
than 0.9, i.e. we exclude pixels that are not completely open to the GRB flux. After 
correcting the arrival time by weighing it with 1/(1 — dead time), we determine from 
the light curve which time intervals have the shape of a FRED distribution. Recall 
from Equation (|8|) that the more photons we take the more we are able to constrain 

K. 

In [37] the average energy difference A(E) = A(E)s — A(E)i was computed for 
each GRB using the energy bands of SWIFT, where A(_E)3 is the average energy of 
the photons with energies between 110 and 300 keV and A(E)i between 20 and 55 
keV (see table [2]) . In our variables the time difference would then be approximately 
given by At = n ■ A(E). 

However, from our analysis we obtain the parameter n directly so we do not 
average over energies in order to get a time difference. Considering only a linear 
approximation to quantum gravitational effects as proposed by Ellis et al 35, 36j, we 
have the relation 



where a and b are coefficients to be fitted. The constant b parameterizes time lags 
in the rest frame of the source caused by unknown internal processes of the GRBs. 
Comparing (J9j) with f4| we find that I(z) is given by 



and the QG parameter a by 



Fitting all GRBs with known redshift detected by INTEGRAL, we find (units s/keV 
are used) 



k = al(z) + 6(1 + z), 



(9) 




(10) 



k = (9.5 ± 3.0) • 10~ 4 • I{z) - (2.8 ± 1.1) • 10" 4 • (1 + z) 



(12) 
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Table 2. Results for the GRBs with known redshifts. Note that in spite of 
the fact that a couple of redshifts have error bars, we choose to take the mean 
value of the redshifts without errors. The reason is that then we don't have 
to introduce arbitrary error bars in order for exact redshifts not to be weighed 
infinitely strongly. I(z) is given by JToJl, K(z) by ||13|I . k is the time lag per energy 
given by the maximization of J5j with cr K its error. 



GRB 


z 


I(z) 


K(z) 


k [s/keV] 




(7 K 


s/keV] 


030227 


1.39 [38] 


1.0 


0.42 


7.8 


10 


4 


4.0 


10" 5 


UOliUO 


nil iQQi 


n in 
U. 1U 


n no 
u.uy 


-1.7 


10- 


4 


n 1 




040106 


0.9 [40] 


0.73 


0.38 


2.7 
4.2 


10" 
10" 


4 
-4 


3.7 
1.5 


10" 4 
10~ 4 


040223 


0.1 gi] 


0.1 


0.09 


2.5 


10" 


3 


0.4 


10 -3 


040812 


0.5 [42] 


0.45 


0.3 


-1.4 
2.6 


10- 

10" 


-3 
-4 


0.1 
0.8 


10" 3 

10- 4 


040827 


0.9 [43] 


0.73 


0.38 


-1.9 


10" 


4 


5.9 


1Q -4 


041218 


0.8 [44] 


0.66 


0.37 


-1.2 
1.7 


10" 
10" 


3 
-3 


0.1 

0.2 


10~ 3 
10~ 3 


050502 


3.8 [S] 


1.69 


0.35 


-2.2 
8.2 


10" 

10- 


3 

-4 


0.4 
1.9 


10 -3 

10- 4 


050714 


0.26 [46] 


0.25 


0.19 


-8.9 


10- 


4 


3.2 


10- 4 


050922 


2.17 g7] 


1.3 


0.41 


-1.4 
7.3 


10" 
10" 


-3 
4 


0.2 
0.2 


10" 3 

10- 4 


060204 


3.1 g8] 


1.55 


0.38 


2.6 
2.3 


10" 
10" 


4 
-4 


6.3 
0.1 


10- 4 
10- 4 



as shown in the left panel of figure 0J Because redshifts are measured without 
using a specific cosmological model, this fit was obtained using data that are model- 
independent. Moreover, a rather questionable energy binning as explained above is 
not needed due to the fact that our analysis method yields directly values for n. 

As can be seen from figure H] (left panel), a single GRB with two bursts can lead 
to very different time lags. For example, GRB040812 with average redshift z — 0.5 has 
two peaks that even differ in the sign: the first one has a negative value k = — 1.4- 1CP 3 
and the second one a positive value n = 2.6 • 10 . This could be explained by the fact 
that different internal processes are at the origin of the two bursts, which implies that 
it may not be sufficient to describe internal time lags with a constant b as in Equation 

However, the physics involved in GRB is still not well understood, thus limiting 
the possibility to model intrinsic effects in other ways than through Equation ©. 

In [36) and [37) a linear fit was obtained by using not z as the independent variable 
but instead a function K(z). Dividing (|9|) by (1 + z), K(z) is given by the non-linear 
function 

K(z) = -±-I{z). (13) 
1 + z 

However, we think that considering k/ (1 + z) as a linear function of K is delusive 
because the new function K(z) is not injective. This function maps certain different 
redshifts z to the same value and has a maximum of i^ max — 0.42 at z = 1.64. For 
example, a redshift of z = 4 has the same value K as a redshift of z = 0.7. Thus the 
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Figure 4. Left panel: plot of function of the redshift z for several GRBs 

detected by INTEGRAL. The nonlinear fit is given by QU, the maximum lies at 
ax — 2.34 with value K max _ ~ 3.6 ■ 10 4 s/keV. Right panel: evolution of the \ 2 



function as a function of M. 
GeV. 



Note that x 2 has a strong minimum around 4 ■ 10 1 



two points for GRB050502 at z — 3.793 are mapped to K = 0.353, which is between 
GRB040812 and GRB040106. Our opinion is that this method is misleading and does 
not give reliable results and should therefore not be used. 

5.2. Likelihood test 

Following [36j . we introduce a likelihood function 

L LH (M)=^exp(-^^) , (14) 

where M is the mass scale, AT the normalization and ^{M) is given by 



2 ,, n v-^ (ftj - b(l + Zj) - a{M)Ii) 
X(M)= ^ . (15) 



all GRBs 

The parameter b reflects the instrinsic time lags and a quantum gravitational effects. 
Thus, b was removed from the linear fit, as can be seen from (|15p . Note that we used 
the raw model that doesn't need an energy binning. 

The value at the minimum of x 2 /d.o.f. is 303/15, which is well above unity. In 
such a case, we may expect a high degree of uncertainty for any fitted parameters. If 
the error bars are underestimated it will lead to underestimated statistical errors for 
the fitted parameters. In such cases, the Particle Data Group [S3] suggests to rescale 
the error bars so that \ 2 ~ d.o.f. by a factor S = [x 2 /d.o.f .J 1 / 2 . Such a rescaling has 
also been proposed in [5SJ [351 13T] 

Figure [4] (right panel) presents the dependence of the rescaled x 2 /d-0.f as a 
function of M . The minimum of this function is found at M ~ 3.8 ■ 10 11 GeV. 
This value also minimizes the likelihood function given by Equation ([5]). 

Following Ellis et al [35] we establish a 95 % confidence-level lower limit on the 
scale M of quantum gravity by solving the equation 

^ £LH(eK = 0.95, (16) 
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where the Planck mass A4 = 10 19 GeV is the reference point fixing the normalization. 
The function L L h is given by Equation (TT4"]) . Solving this equation for M gives the 
lower limit of quantum gravity at a 95 % level of confidence at 

M > 3.2 • 10 11 GeV. (17) 
6. Conclusions 

In this work, we first described a method that is able to analyze unbinned data of GRBs 
detected by INTEGRAL. We introduced a maximum likelihood function following a 
Fast Raise and Exponential Decay behavior with a parameter describing time lags 
of photons for different energies. In order to know which minimum time lags are 
measurable with INTEGRAL, we performed Monte Carlo simulations and varied the 
total photon number. 

We had 11 GRBs with known redshift at our disposal and were able to get 17 
measurements of time lags. We used these measurements to fit a nonlinear relation 
depending on the redshift. This relation has a term that describes possible quantum 
gravitational effects and one that accounts for intrinsic time lags of the GRB. By 
using a likelihood function we made a \ 2 analysis of the data and showed that there 
is a strong minimum of % 2 around 4 • 10 11 GeV, which apparently would disfavor a 
quantum gravitational scale around the Planck mass. However, as shown by our Monte 
Carlo simulations in section [4] it is obvious that it is impossible to obtain the required 
sensitivity with the presently available statistics of GRB data, especially when only 
11 GRBs are at disposal. Correcting for intrisinc time lags [55J (56J, [57] dramatically 
increases this lower bound to 1.5 • 10 14 GeV, but this method stands on shaky ground. 

A better precision in time could be achieved by constructing satellites with a much 
larger collecting surface. However, as shown in section [4] even with the unrealistic 
photon number of 3 • 10 5 for a single burst the time resolution is still two orders of 
magnitude too low. The other solution is to increase the photon energy, as can be 
seen from Equation (g]) . A more complete strategy how to reach the strongest possible 
bounds was discussed in [5S1 [SHj ■ GLAST will be able to increase the time resolution 
by several orders of magnitude as it will be able to detect photons up to energies of 20 
MeV. As the time lags are linearly dependent on the energy difference GLAST should 
be able to improve the time difference by a couple of orders of magnitude, which may 
be even larger than the expected time difference caused by QG effects at the Planck 
scale. 
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